Spatial transcriptome analysis of long non-coding RNAs reveals tissue specificity and functional roles in cancer / 浙江大学学报（英文版）（B辑：生物医学和生物技术）

Long non-coding RNAs (lncRNAs) play a significant role in maintaining tissue morphology and functions, and their precise regulatory effectiveness is closely related to expression patterns. However, the spatial expression patterns of lncRNAs in humans are poorly characterized. Here, we constructed five comprehensive transcriptomic atlases of human lncRNAs covering thousands of major tissue samples in normal and disease states. The lncRNA transcriptomes exhibited high consistency within the same tissues across resources, and even higher complexity in specialized tissues. Tissue-elevated (TE) lncRNAs were identified in each resource and robust TE lncRNAs were refined by integrative analysis. We detected 1 to 4684 robust TE lncRNAs across tissues; the highest number was in testis tissue, followed by brain tissue. Functional analyses of TE lncRNAs indicated important roles in corresponding tissue-related pathways. Moreover, we found that the expression features of robust TE lncRNAs made them be effective biomarkers to distinguish tissues; TE lncRNAs also tended to be associated with cancer, and exhibited differential expression or were correlated with patient survival. In summary, spatial classification of lncRNAs is the starting point for elucidating the function of lncRNAs in both maintenance of tissue morphology and progress of tissue-constricted diseases.

Amino acid metabolism characteristics of Banxia Baizhu Tianma Decoction in realizing drug withdrawal based on transcriptomic analysis / 中国中药杂志

This study aimed to demonstrate the effect of Banxia Baizhu Tianma Decoction(BBTD) on realizing withdrawal of anti-epileptic drugs and explore the relationship between BBTD and the amino acid metabolism by transcriptomic analysis in the rat model of epilepsy induced by lithium chloride-pilocarpine. The rats with epilepsy were divided into a control group(Ctrl), an epilepsy group(Ep), a BBTD & antiepileptic drug integrative group(BADIG), and an antiepileptic drug withdrawal group(ADWG). The Ctrl and Ep were given ultrapure water by gavage for 12 weeks. The BADIG was given BBTD extract and carbamazepine solution by gavage for 12 weeks. The ADWG was given carbamazepine solution and BBTD extract by gavage for the former 6 weeks, and then only given BBTD extract for the latter 6 weeks. The therapeutic effect was evaluated by behavioral observation, electroencephalogram(EEG), and hippocampal neuronal morphological changes. High-throughput sequencing was used to obtain amino acid metabolism-related differen-tial genes in the hippocampus, and the mRNA expression in the hippocampus of each group was verified by real-time quantitative polymerase chain reaction(RT-qPCR). The hub genes were screened out through protein-protein interaction(PPI) network, and Gene Ontology(GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment analysis were performed. Two ceRNA networks, namely circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA, were constructed for ADWG vs BADIG. The experimental results showed that compared with those in Ep, rats in ADWG were significantly improved in the behavioral observation, EEG, and hippocampal neuronal impairment. Thirty-four amino acid metabolism-related differential genes were obtained by transcriptomic analysis, and the sequencing results were confirmed by RT-qPCR. Eight hub genes were obtained through PPI network, involving several biological processes, molecular functions, and signal pathways related to amino acid metabolism. Finally, the circRNA-miRNA-mRNA ternary transcription network of 17 circRNA, 5 miRNA, and 2 mRNA, and a lncRNA-miRNA-mRNA ternary network of 10 lncRNA, 5 miRNA, and 2 mRNA were constructed in ADWG vs BADIG. In conclusion, BBTD can effectively achieve the withdrawal of antiepileptic drugs, which may be related to the transcriptomic regulation of amino acid metabolism.

Comparative analysis of differentially expressed genes for biosynthesis of active ingredients in fruits of different cultivars of Lycium barbarum L. based on transcriptome sequencing / 生物工程学报

To explore the differentially expressed genes (DEGs) related to biosynthesis of active ingredients in wolfberry fruits of different varieties of Lycium barbarum L. and reveal the molecular mechanism of the differences of active ingredients, we utilized Illumina NovaSeq 6000 high-throughput sequencing technology to conduct transcriptome sequencing on the fruits of 'Ningqi No.1' and 'Ningqi No.7' during the green fruit stage, color turning stage and maturity stage. Subsequently, we compared the profiles of related gene expression in the fruits of the two varieties at different development stages. The results showed that a total of 811 818 178 clean reads were obtained, resulting in 121.76 Gb of valid data. There were 2 827, 2 552 and 2 311 DEGs obtained during the green fruit stage, color turning stage and maturity stage of 'Ningqi No. 1' and 'Ningqi No. 7', respectively, among which 2 153, 2 050 and 1 825 genes were annotated in six databases, including gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG) and clusters of orthologous groups of proteins (KOG). In GO database, 1 307, 865 and 624 DEGs of green fruit stage, color turning stage and maturity stage were found to be enriched in biological processes, cell components and molecular functions, respectively. In the KEGG database, the DEGs at three developmental stages were mainly concentrated in metabolic pathways, biosynthesis of secondary metabolites and plant-pathogen interaction. In KOG database, 1 775, 1 751 and 1 541 DEGs were annotated at three developmental stages, respectively. Searching the annotated genes against the PubMed database revealed 18, 26 and 24 DEGs related to the synthesis of active ingredients were mined at the green fruit stage, color turning stage and maturity stage, respectively. These genes are involved in carotenoid, flavonoid, terpenoid, alkaloid, vitamin metabolic pathways, etc. Seven DEGs were verified by RT-qPCR, which showed consistent results with transcriptome sequencing. This study provides preliminary evidences for the differences in the content of active ingredients in different Lycium barbarum L. varieties from the transcriptional level. These evidences may facilitate further exploring the key genes for active ingredients biosynthesis in Lycium barbarum L. and analyzing their expression regulation mechanism.

Identification and expression pattern analysis of RcACA gene family in castor under abiotic stresses / 生物工程学报

Auto-inhibited Ca2+-ATPase (ACA) is one of the Ca2+-ATPase subfamilies that plays an important role in maintaining Ca2+ concentration balance in plant cells. To explore the function and gene expression pattern of the RcACA gene family in castor, bioinformatics analysis was used to identify the members of the RcACA gene family in castor. The basic physical and chemical properties, subcellular location, protein secondary and tertiary structure, conserved domain, conserved motif, gene structure, chromosome location and collinear relationship, as well as the evolutionary characteristics and promoter cis-acting elements were predicted and analyzed. The expression pattern of the RcACA gene under abiotic stress was analyzed by expression (fragments per kilobase of exon model per million mapped fragments, FPKM) in castor transcriptome data. The results showed that 8 RcACA gene family members were identified in castor, acidic proteins located in the plasma membrane. In the secondary structure of all proteins, the α-helix and random coil is more; the RcACA genes were clustered into three categories, and the design of the genes in the same category was similar to the conserved motif. Both of them had four typical domains, RcACA3-RcACA8 had a Ca2+-ATPase N-terminal autoinhibitory domain. The RcACA gene is mostly located on the long arm of the chromosome and has 2 pairs of collinear relationships. There are more light response elements but fewer hormone-induced elements located upstream of the RcACA coding region. Interspecific clustering showed that the evolution of ACA genes among species was conservative. Tissue expression pattern analysis showed that RcACA genes showed apparent tissue expression specificity, and most of the genes showed the highest expression level in male flowers. Expression analysis under abiotic stress showed that RcACA2-RcACA8 were up-regulated under high salt and drought stress, and RcACA1 was up-regulated at 0-24 h under low-temperature stress, indicating that RcACA genes positively responded to abiotic stresses. The above results provide a theoretical basis for exploring the role of the RcACA gene in castor growth, development and stress response.

The role of B3GNT3 as an oncogene in the growth, invasion and migration of esophageal cancer cells

Purpose: To investigate the role and mechanism of ß1,3-N-acetylglucosaminyltransferase-3 gene (B3GNT3) in esophageal cancer (ESCA). Methods: The starBase database was used to evaluate the expression of B3GNT3. B3GNT3 function was measured using KYSE-30 and KYSE-410 cells of esophageal squamous cell carcinoma (ESCC) cell lines. The mRNA levels were detected by quantitative real-time polymerase chain reaction (qRT-PCR). Cell counting kit-8, clone formation assay and transwell assay were used to detect the changes of proliferation, invasion and migration. Results: B3GNT3 expression was higher in ESCA tissues than in normal tissues. The overall survival rate of ESCA patients with high B3GNT3 expression was lower than that of ESCA patients with low B3GNT3 expression. In vitro functional experiments showed that the proliferation ability, migration and invasion ability of KYSE-30 and KYSE-410 cells with B3GNT3 interference were lower than those of the control, and the overexpression of B3GNT3 had the opposite effect. After silencing B3GNT3 expression in ESCC cell lines, the growth of both cell lines was inhibited and the invasiveness was decreased. Knockdown of B3GNT3 reduced the growth rate and Ki-67 expression level. Conclusion: B3GNT3, as an oncogene, may promote the growth, invasion and migration of ESCC cell.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L) response to repeated drought stresses

BACKGROUND: Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions. RESULTS: We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes. CONCLUSION: Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

Single-cell RNA sequencing reveals the transcriptomic landscape of kidneys in patients with ischemic acute kidney injury / 中华医学杂志(英文版)

BACKGROUND@#Ischemic acute kidney injury (AKI) is a common syndrome associated with considerable mortality and healthcare costs. Up to now, the underlying pathogenesis of ischemic AKI remains incompletely understood, and specific strategies for early diagnosis and treatment of ischemic AKI are still lacking. Here, this study aimed to define the transcriptomic landscape of AKI patients through single-cell RNA sequencing (scRNA-seq) analysis in kidneys.@*METHODS@#In this study, scRNA-seq technology was applied to kidneys from two ischemic AKI patients, and three human public scRNA-seq datasets were collected as controls. Differentially expressed genes (DEGs) and cell clusters of kidneys were determined. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, as well as the ligand-receptor interaction between cells, were performed. We also validated several DEGs expression in kidneys from human ischemic AKI and ischemia/reperfusion (I/R) injury induced AKI mice through immunohistochemistry staining.@*RESULTS@#15 distinct cell clusters were determined in kidney from subjects of ischemic AKI and control. The injured proximal tubules (PT) displayed a proapoptotic and proinflammatory phenotype. PT cells of ischemic AKI had up-regulation of novel pro-apoptotic genes including USP47 , RASSF4 , EBAG9 , IER3 , SASH1 , SEPTIN7 , and NUB1 , which have not been reported in ischemic AKI previously. Several hub genes were validated in kidneys from human AKI and renal I/R injury mice, respectively. Furthermore, PT highly expressed DEGs enriched in endoplasmic reticulum stress, autophagy, and retinoic acid-inducible gene I (RIG-I) signaling. DEGs overexpressed in other tubular cells were primarily enriched in nucleotide-binding and oligomerization domain (NOD)-like receptor signaling, estrogen signaling, interleukin (IL)-12 signaling, and IL-17 signaling. Overexpressed genes in kidney-resident immune cells including macrophages, natural killer T (NKT) cells, monocytes, and dendritic cells were associated with leukocyte activation, chemotaxis, cell adhesion, and complement activation. In addition, the ligand-receptor interactions analysis revealed prominent communications between macrophages and monocytes with other cells in the process of ischemic AKI.@*CONCLUSION@#Together, this study reveals distinct cell-specific transcriptomic atlas of kidney in ischemic AKI patients, altered signaling pathways, and potential cell-cell crosstalk in the development of AKI. These data reveal new insights into the pathogenesis and potential therapeutic strategies in ischemic AKI.

EHHADH is a key gene in fatty acid metabolism pathways in hepatocellular carcinoma: a transcriptomic analysis / 南方医科大学学报

OBJECTIVE@#To explore the driving gene of hepatocellular carcinoma (HCC) occurrence and progression and its potential as new therapeutic target of HCC.@*METHODS@#The transcriptome and genomic data of 858 HCC tissues and 493 adjacent tissues were obtained from TCGA, GEO, and ICGC databases. Gene Set Enrichment Analysis (GSEA) identified EHHADH (encoding enoyl-CoA hydratase/L-3-hydroxyacyl-CoA dehydrogenase) as the hub gene in the significantly enriched differential pathways in HCC. The downregulation of EHHADH expression at the transcriptome level was found to correlate with TP53 mutation based on analysis of the TCGA- HCC dataset, and the mechanism by which TP53 mutation caused EHHADH downregulation was explored through correlation analysis. Analysis of the data from the Metascape database suggested that EHHADH was strongly correlated with the ferroptosis signaling pathway in HCC progression, and to verify this result, immunohistochemical staining was used to examine EHHADH expression in 30 HCC tissues and paired adjacent tissues.@*RESULTS@#All the 3 HCC datasets showed signficnatly lowered EHHADH expression in HCC tissues as compared with the adjacent tissues (P < 0.05) with a close correlation with the degree of hepatocyte de-differentiation (P < 0.01). The somatic landscape of HCC cohort in TCGA dataset showed that HCC patients had the highest genomic TP53 mutation rate. The transcriptomic level of PPARGC1A, the upstream gene of EHHADH, was significantly downregulated in HCC patients with TP53 mutation as compared with those without the mutation (P < 0.05), and was significantly correlated with EHHADH expression level. GO and KEGG enrichment analyses showed that EHHADH expression was significantly correlated with abnormal fatty acid metabolism in HCC. The immunohistochemical results showd that the expression level of EHHADH in HCC tissues was down-regulated, and its expression level was related to the degree of hepatocytes de-differentiation and the process of ferroptosis.@*CONCLUSION@#TP53 mutations may induce abnormal expression of PPARGC1A to cause downregulation of EHHADH expression in HCC. The low expression of EHHADH is closely associated with aggravation of de-differentiation and ferroptosis escape in HCC tissues, suggesting the potential of EHHADH as a therapeutic target for HCC.

Identification of core pathogenic genes and pathways in elderly osteoporosis based on bioinformatics analysis / 中华预防医学杂志

Objective: Using bioinformatics methods to analyze the core pathogenic genes and related pathways in elderly osteoporosis. Methods: From November 2020 and August 2021, eight elderly osteoporosis patients who received treatment and five healthy participants who underwent physical examinations in Beijing Jishuitan Hospital were selected as subjects. The expression level of RNA in the peripheral blood of eight elderly osteoporosis patients and five healthy participants was collected for high-throughput transcriptome sequencing and analysis. The gene ontology (GO) analysis Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed for the differentially expressed genes (DEGs). The protein-protein interaction (PPI) network was constructed using the STRING website and Cytoscape software, and the most significant modules and hub genes were screened out. Results: Among the eight elderly osteoporosis patients, there were seven females and one male, with an average age of 72.4 years (SD=4.2). Among the five healthy participants, there were four females and one male, with an average age of 68.2 years (SD=5.7). A total of 1 635 DEGs (847 up-regulated and 788 down-regulated) were identified. GO analysis revealed that the molecular functions of DEGs were mainly enriched in structural constituents of the ribosome, protein dimerization activity, and cellular components were mainly enriched in the nucleosome, DNA packaging complex, cytosolic part, protein-DNA complex and the cytosolic ribosome. KEGG pathway analysis showed that DEGs were mainly enriched in systemic lupus erythematosus and ribosome. Gene UBA52, UBB, RPS27A, RPS15, RPS12, RPL13A, RPL23A, RPL10A, RPS25 and RPS6 were selected and seven of them could encode ribosome proteins. Conclusion: The pathogenesis of elderly osteoporosis may be associated with ribosome-related genes and pathways.

Developmental effects of TCIPP and TnBP on zebrafish (Danio rerio) embryos / 中华预防医学杂志

Objective: To investigate the toxicity of tris (2-chloropropyl) phosphate (TCIPP) and tributyl phosphate (TnBP) on the growth and development of zebrafish embryos, as well as to explore the underlying mechanisms at the transcriptional level. Methods: With zebrafish as a model, two hpf zebrafish embryos were exposed to TCIPP and TnBP (0.1, 1, 10, 100, 500, and 1 000 μmol/L) using the semi-static method, and their rates of lethality and hatchability were determined. The transcriptome changes of 120 hpf juvenile zebrafish exposed to environmentally relevant concentrations of 0.1 and 1 μmol/L were measured. Results: The 50% lethal concentrations (LC50) of TCIPP and TnBP for zebrafish embryos were 155.30 and 27.62 μmol/L (96 hpf), 156.5 and 26.05 μmol/L (120 hpf), respectively. The 72 hpf hatching rates of TCIPP (100 μmol/L) and TnBP (10 μmol/L) were (23.33±7.72)% and (91.67±2.97)%, which were significantly decreased compared with the control group (P<0.05). Transcriptome analysis showed that TnBP had more differential genes (DEGs) than TCIPP, with a dose-response relationship. These DEGs were enriched in 32 pathways in total, including those involved in oxidative stress, energy metabolism, lipid metabolism, and nuclear receptor-related pathways, using the IPA pathway analysis. Among them, three enriched pathways overlapped between TCIPP and TnBP, including TR/RXR activation and CAR/RXR activation. Additionally, DEGs were also mapped onto pathways of LXR/RXR activation and oxidative stress for TnBP exposure only. Conclusion: Both TCIPP and TnBP have growth and developmental toxicities in zebrafish embryos, with distinct biomolecular mechanisms, and TnBP has a stronger effect than TCIPP.

Mechanism of artesunate on bone destruction in experimental rheumatoid arthritis based on transcriptomics and network pharmacology / 中国中药杂志

The present study investigated the mechanism of artesunate in the treatment of bone destruction in experimental rheumatoid arthritis(RA) based on transcriptomics and network pharmacology. The transcriptome sequencing data of artesunate in the inhibition of osteoclast differentiation were analyzed to obtain differentially expressed genes(DEGs). GraphPad Prism 8 software was used to plot volcano maps and heat maps were plotted through the website of bioinformatics. GeneCards and OMIM were used to collect information on key targets of bone destruction in RA. The DEGs of artesunate in inhibiting osteoclast differentiation and key target genes of bone destruction in RA were intersected by the Venny 2.1.0 platform, and the intersection target genes were analyzed by Gene Ontology(GO)/Kyoto Encyclopedia of Genes and Genomes(KEGG) enrichment. Finally, the receptor activator of nuclear factor-κB(RANKL)-induced osteoclast differentiation model and collagen-induced arthritis(CIA) model were established. Quantitative real time polymerase chain reaction(q-PCR), immunofluorescence, and immunohistochemistry were used to verify the pharmacological effect and molecular mechanism of artesunate in the treatment of bone destruction in RA. In this study, the RANKL-induced osteoclast differentiation model in vitro was established and intervened with artesunate, and transcriptome sequencing data were analyzed to obtain 744 DEGs of artesunate in inhibiting osteoclast differentiation. A total of 1 291 major target genes of bone destruction in RA were obtained from GeneCards and OMIM. The target genes of artesunate in inhibiting osteoclast differentiation and the target genes of bone destruction in RA were intersected to obtain 61 target genes of artesunate against bone destruction in RA. The intersected target genes were analyzed by GO/KEGG enrichment. According to the results previously reported, the cytokine-cytokine receptor interaction signaling pathway was selected for experimental verification. Artesunate intervention in the RANKL-induced osteoclast differentiation model showed that artesunate inhibited CC chemokine receptor 3(CCR3), CC chemokine receptor 1(CCR1) and leukemia inhibitory factor(LIF) mRNA expression in osteoclasts in a dose-dependent manner compared with the RANKL-induced group. Meanwhile, the results of immunofluorescence and immunohistochemistry showed that artesunate could dose-dependently reduce the expression of CCR3 in osteoclasts and joint tissues of the CIA rat model in vitro. This study indicated that artesunate regulated the CCR3 in the cytokine-cytokine receptor interaction signaling pathway in the treatment of bone destruction in RA and provided a new target gene for the treatment of bone destruction in RA.

Mechanism of Huangjing Qianshi Decoction in treatment of prediabetic mice based on transcriptome sequencing / 中国中药杂志

Based on transcriptome sequencing technology, the mouse model of prediabetes treated with Huangjing Qianshi Decoction was sequenced to explore the possible mechanism of treating prediabetes. First of all, transcriptome sequencing was performed on the normal BKS-DB mouse group, the prediabetic model group, and the Huangjing Qianshi Decoction treatment group(treatment group) to obtain differentially expressed genes in the skeletal muscle samples of mice. The serum biochemical indexes were detected in each group to screen out the core genes of Huangjing Qianshi Decoction in prediabetes. Gene Ontology(GO) database and Kyoto Encyclopedia of Genes and Genomes(KEGG) database were used to conduct signaling pathway enrichment analysis of differentially expressed genes, and real-time quantitative polymerase chain reaction(RT-qPCR) was used to verify them. The results showed that the levels of fasting blood glucose(FBG), fasting insulin(FINS), insulin resistance index(HOMA-IR), total cholesterol(TC), triglycerides(TG), and low-density lipoprotein cholesterol(LDL-C) in the mouse model were significantly decreased after treatment with Huangjing Qianshi Decoction. In the results of differential gene screening, there were 1 666 differentially expressed genes in the model group as compared with the normal group, and there were 971 differentially expressed genes in the treatment group as compared with the model group. Among them, interleukin-6(IL-6) and NR3C2 genes, which were closely related to the regulation of insulin resis-tance function, were significantly up-regulated between the model group and the normal group, and vascular endothelial growth factor A(VEGFA) genes were significantly down-regulated between the model group and the normal group. However, the expression results of IL-6, NR3C2, and VEGFA genes were adverse between the treatment group and the model group. GO functional enrichment analysis found that the biological process annotation mainly focused on cell synthesis, cycle, and metabolism; cell component annotation mainly focused on organelles and internal components; and molecular function annotation mainly focused on binding molecular functions. KEGG pathway enrichment analysis found that it involved the protein tyrosine kinase 6(PTK6) pathway, CD28-dependent phosphoinositide 3-kinase/protein kinase B(PI3K/AKT) pathway, p53 pathway, etc. Therefore, Huangjing Qianshi Decoction can improve the state of prediabetes, and the mechanism may be related to cell cycle and apoptosis, PI3K/AKT pathway, p53 pathway, and other biological pathways regulated by IL-6, NR3C2, and VEGFA.

Analysis of WRKY transcription factor family based on full-length transcriptome sequencing in Polygonatum cyrtonema / 中国中药杂志

WRKY transcription factor family plays an important role in plant growth and development, secondary metabolite synthesis, and biotic and abiotic stress responses. The present study performed full-length transcriptome sequencing of Polygonatum cyrtonema by virtue of the PacBio SMRT high-throughput platform, identified the WRKY family by bioinformatics methods, and analyzed the physicochemical properties, subcellular localization, phylogeny, and conserved motifs. The results showed that 30.69 Gb nucleotide bases and 89 564 transcripts were obtained after redundancy removal. These transcripts had a mean length of 2 060 bp and an N50 value of 3 156 bp. Based on the full-length transcriptome sequencing data, 64 candidate proteins were selected from the WRKY transcription factor family, with the protein size of 92-1 027 aa, the relative molecular mass of 10 377.85-115 779.48 kDa, and the isoelectric point of 4.49-9.84. These WRKY family members were mostly located in the nucleus and belonged to the hydrophobic proteins. According to the phylogenetic analysis of WRKY family in P. cyrtonema and Arabidopsis thaliana, all WRKY family members were clustered into seven subfamilies and WRKY proteins from P. cyrtonema were distributed in different numbers in these seven subgroups. Expression pattern analysis confirmed that 40 WRKY family members had distinct expression patterns in the rhizomes of 1-and 3-year-old P. cyrtonema. Except for PcWRKY39, the expression of 39 WRKY family members was down-regulated in 3-year-old samples. In conclusion, this study provides abundant reference data for genetic research on P. cyrtonema and lays a foundation for the in-depth investigation of the biological functions of the WRKY family.

Preliminary Results of Single-cell Transcriptome Sequencing in Renal Arterial Lesions of Takayasu Arteritis / 中国医学科学院学报

Objective To explore the preliminary application of single-cell RNA sequencing (scRNA-seq) in the renal arterial lesions in Takayasu arteritis (TA) patients. Methods This study included 2 TA patients with renal artery stenosis treated by bypass surgery in the Department of Vascular Surgery,Beijing Hospital.The obtained 2 renal artery samples were digested with two different protocols (GEXSCOPE kit and self-made digestion liquid) before scRNA-seq and bioinformatics analysis. Results A total of 2920 cells were obtained for further analysis.After unbiased cluster analysis,2 endothelial cell subsets,2 smooth muscle cell subsets,1 fibroblast subset,2 mononuclear macrophage subsets,1 T cell subset,and 1 undefined cell subset were identified.Among them,the two subsets of smooth muscle cells were contractile and secretory,respectively.The results of scRNA-seq indicated that enzymatic hydrolysis with GEXSCOPE kit produced a large number of endothelial cells (57.46%) and a small number of immune cells (13.21%).However,immune cells (34.64%) were dominant in the cells obtained by enzymatic hydrolysis with self-made digestive liquid. Conclusion scRNA-seq can be employed to explore the cellular heterogeneity of diseased vessels in TA patients.Different enzymatic digestion protocols may impact the proportion of different cells.

Screening of housekeeping genes in Gelsemium elegans and expression patterns of genes involved in its alkaloid biosynthesis / 生物工程学报

Gelsemium elegans is a traditional Chinese herb of medicinal importance, with indole terpene alkaloids as its main active components. To study the expression of the most suitable housekeeping reference genes in G. elegans, the root bark, stem segments, leaves and inflorescences of four different parts of G. elegans were used as materials in this study. The expression stability of 10 candidate housekeeping reference genes (18S, GAPDH, Actin, TUA, TUB, SAND, EF-1α, UBC, UBQ, and cdc25) was assessed through real-time fluorescence quantitative PCR, GeNorm, NormFinder, BestKeeper, ΔCT, and RefFinder. The results showed that EF-1α was stably expressed in all four parts of G. elegans and was the most suitable housekeeping gene. Based on the coexpression pattern of genome, full-length transcriptome and metabolome, the key candidate targets of 18 related genes (AS, AnPRT, PRAI, IGPS, TSA, TSB, TDC, GES, G8H, 8-HGO, IS, 7-DLS, 7-DLGT, 7-DLH, LAMT, SLS, STR, and SGD) involved in the Gelsemium alkaloid biosynthesis were obtained. The expression of 18 related enzyme genes were analyzed by qRT-PCR using the housekeeping gene EF-1α as a reference. The results showed that these genes' expression and gelsenicine content trends were correlated and were likely to be involved in the biosynthesis of the Gelsemium alkaloid, gelsenicine.

Single-nucleus transcriptomics reveals a gatekeeper role for FOXP1 in primate cardiac aging

Aging poses a major risk factor for cardiovascular diseases, the leading cause of death in the aged population. However, the cell type-specific changes underlying cardiac aging are far from being clear. Here, we performed single-nucleus RNA-sequencing analysis of left ventricles from young and aged cynomolgus monkeys to define cell composition changes and transcriptomic alterations across different cell types associated with age. We found that aged cardiomyocytes underwent a dramatic loss in cell numbers and profound fluctuations in transcriptional profiles. Via transcription regulatory network analysis, we identified FOXP1, a core transcription factor in organ development, as a key downregulated factor in aged cardiomyocytes, concomitant with the dysregulation of FOXP1 target genes associated with heart function and cardiac diseases. Consistently, the deficiency of FOXP1 led to hypertrophic and senescent phenotypes in human embryonic stem cell-derived cardiomyocytes. Altogether, our findings depict the cellular and molecular landscape of ventricular aging at the single-cell resolution, and identify drivers for primate cardiac aging and potential targets for intervention against cardiac aging and associated diseases.

Single-cell transcriptomic atlas of mouse cochlear aging

Progressive functional deterioration in the cochlea is associated with age-related hearing loss (ARHL). However, the cellular and molecular basis underlying cochlear aging remains largely unknown. Here, we established a dynamic single-cell transcriptomic landscape of mouse cochlear aging, in which we characterized aging-associated transcriptomic changes in 27 different cochlear cell types across five different time points. Overall, our analysis pinpoints loss of proteostasis and elevated apoptosis as the hallmark features of cochlear aging, highlights unexpected age-related transcriptional fluctuations in intermediate cells localized in the stria vascularis (SV) and demonstrates that upregulation of endoplasmic reticulum (ER) chaperon protein HSP90AA1 mitigates ER stress-induced damages associated with aging. Our work suggests that targeting unfolded protein response pathways may help alleviate aging-related SV atrophy and hence delay the progression of ARHL.

Multiplexed single-cell transcriptome analysis reveals molecular characteristics of monkey pluripotent stem cell lines / 浙江大学学报（英文版）（B辑：生物医学和生物技术）

Efforts have been made to establish various human pluripotent stem cell lines. However, such methods have not yet been duplicated in non-human primate cells. Here, we introduce a multiplexed single-cell sequencing technique to profile the molecular features of monkey pluripotent stem cells in published culture conditions. The results demonstrate suboptimized maintenance of pluripotency and show that the selected signaling pathways for resetting human stem cells can also be interpreted for establishing monkey cell lines. Overall, this work legitimates the translation of novel human cell line culture conditions to monkey cells and provides guidance for exploring chemical cocktails for monkey stem cell line derivation.

Genome-wide DNA methylation and transcriptome expression profiles of peripheral blood mononuclear cells in patients with systemic sclerosis with interstitial lung disease / 中南大学学报(医学版)

OBJECTIVES@#This study aims to investigate the genome-wide DNA methylation and transcriptome expression profiles of peripheral blood mononuclear cells (PBMCs) in patients with systemic sclerosis (SSc) with interstitial lung disease (ILD), and to analyze the effects of DNA methylation on Wnt/β-catenin and chemokine signaling pathways.@*METHODS@#PBMCs were collected from 19 patients with SSc (SSc group) and 18 healthy persons (control group). Among SSc patients, there were 10 patients with ILD (SSc with ILD subgroup) and 9 patients without ILD (SSc without ILD subgroup). The genome-wide DNA methylation and gene expression level were analyzed by using Illumina 450K methylation chip and Illumina HT-12 v4.0 gene expression profiling chip. The effect of DNA methylation on Wnt/β-catenin and chemokine signal pathways was investigated.@*RESULTS@#Genome-wide DNA methylation analysis identified 71 hypermethylated CpG sites and 98 hypomethylated CpG sites in the SSc with ILD subgroup compared with the SSc without ILD subgroup. Transcriptome analysis distinguished 164 upregulated genes and 191 downregulated genes in the SSc with ILD subgroup as compared with the SSc without ILD subgroup. In PBMCs of the SSc group, 35 genes in Wnt/β-catenin signaling pathway were hypomethylated, while frizzled-1 (FZD1), mitogen-activated protein kinase 9 (MAPK9), mothers against DPP homolog 2 (SMAD2), transcription factor 7-like 2 (TCF7L2), and wingless-type MMTV integration site family, member 5B (WNT5B) mRNA expressions were upregulated as compared with the control group (all P<0.05). Compared with the SSc without ILD subgroup, the mRNA expressions of dickkopf homolog 2 (DKK2), FZD1, MAPK9 were upregulated in the SSc with ILD subgroup, but the differences were not statistically significant (all P>0.05). In PBMCs of the SSc group, 38 genes in chemokine signaling pathway were hypomethylated, while β-arrestin 1 (ARRB1), C-X-C motif chemokine ligand 10 (CXCL10), C-X-C motif chemokine ligand 16 (CXCL16), FGR, and neutrophil cytosolic factor 1C (NCF1C) mRNA expressions were upregulated as compared with the control group (all P<0.05). Compared with the SSc without ILD subgroup, the mRNA expressions of ARRB1, CXCL10, CXCL16 were upregulated in the SSc with ILD subgroup, but the differences were not statistically significant (all P>0.05).@*CONCLUSIONS@#There are differences in DNA methylation and transcriptome profiles between SSc with ILD and SSc without ILD. The expression levels of multiple genes in Wnt/β- catenin and chemokine signaling pathways are upregulated, which might be associatea with the pathogenesis of SSc.

Shexiang Tongxin Dropping Pill Allieviates Heart Failure via Extracellula Matrix-Receptor Interaction Pathways Based on RNA-Seq Transcriptomics and Experimental Studies / 中国结合医学杂志

OBJECTIVE@#To investigate the protective mechanisms of Chinese medicine Shexiang Tongxin Dropping Pills (STDP) on heart failure (HF).@*METHODS@#Isoproterenol (ISO)-induced HF rat model and angiotensin II (Ang II)-induced neonatal rat cardiac fibroblast (CFs) model were used in the present study. HF rats were treated with and without STDP (3 g/kg). RNA-seq was performed to identify differentially expressed genes (DEGs). Cardiac function was evaluated by echocardiography. Hematoxylin and eosin and Masson's stainings were taken to assess cardiac fibrosis. The levels of collagen I (Col I) and collagen III (Col III) were detected by immunohistochemical staining. CCK8 kit and transwell assay were implemented to test the CFs' proliferative and migratory activity, respectively. The protein expressions of α-smooth muscle actin (α-SMA), matrix metalloproteinase-2 (MMP-2), MMP-9, Col I, and Col III were detected by Western blotting.@*RESULTS@#The results of RNA-seq analysis showed that STDP exerted its pharmacological effects on HF via multiple signaling pathways, such as the extracellular matrix (ECM)-receptor interaction, cell cycle, and B cell receptor interaction. Results from in vivo experiments demonstrated that STDP treatment reversed declines in cardiac function, inhibiting myocardial fibrosis, and reversing increases in Col I and Col III expression levels in the hearts of HF rats. Moreover, STDP (6, 9 mg/mL) inhibited the proliferation and migration of CFs exposed to Ang II in vitro (P<0.05). The activation of collagen synthesis and myofibroblast generation were markedly suppressed by STDP, also the synthesis of MMP-2 and MMP-9, as well as ECM components Col I, Col III, and α-SMA were decreased in Ang II-induced neonatal rats' CFs.@*CONCLUSIONS@#STDP had anti-fibrotic effects in HF, which might be caused by the modulation of ECM-receptor interaction pathways. Through the management of cardiac fibrosis, STDP may be a compelling candidate for improving prognosis of HF.